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ABSTRACT 

In this, the first of a series of papers, we investigate the variation 
along a magnetic field line of the energy and pitch angle distribution of high 
energy electrons injected into a cold hydrogen plasma containing either an 
open or closed magnetic field structure. The problem is formulated as a time- 
independent Fokker-Planck Equation for the electron number distribution as a 
function of the electron energy s electron pitch angle and the structure of the 
global magnetic field. 

We present a simple analytic solution valid in the small pitch angle regime 
and for a slowly varying magnetic field. For the more general situation we 
use a numerical code for solving the Fokker-Planck Equation, and we find that 
the analytic expression agrees v^/ell with the numerical results to values of 
the pitch angle much larger than expocted. For most practical applications 
one may confidently use the analytic expression instead of having to resort to 
lengthy numerical computations. 

Those results may be useful for a variety of astrophysieal applications. 

Our primary concern is to use these results to study the non -thermal models of 
the impulsive phase of solar flares. In subsequent papers we shall use the 
results of this paper to calculate the expected x-ray and microwave radiation 
from a flaring loop and compare those results with the high resolution data 
currently becoming available. 

Subject He adings : Atomic processes, Sun: Flares, Sun: X-rays, Sun:Radio Radiation 
*Also Department of Applied Physics 
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I. INTRODUCTION 

Hard x-rays and microwave radiation observed during the Impulsive phase of a 
solar flare provide us with a means for Inferring the characteristics of the 
distribution of energetic electrons within the flaring plasma. Although It Is 
now widely accepted that the hard x-rays and microwaves are the result of 
coulomb bremsstrahlung and gyrosynchrotron mechanisms respectively, there Is 
still no consensus as to the form of the electron distribution function within 
the flaring loop nor as to the relative Importance of the various candidate 
processes which may affect or determine the characteristics of this distribution. 

Various models describing the operative physical processes have been 
proposed, broadly classified as thermal or non-thermal ntodels (see Emslie and 
Rust 1979 and references therein). The non-thennal models may be subdivided 
into thick or thin target beamed electron models and trapped electron models. 

The prime differences between these models have often been taken to lie jwith 
the spectral distribution of accelerated electrons and the characteristics of 
the background plasma. (For a discussion of these models the interested reader 
is referred to Sturrock 1980, p. 21 5, and to Brown and Smith 1980). 

Principally two types of mechanism have been suggested by which a situation 
giving rise to the emission of thermal hard x-ray bremsstrahlung from a flare 
may be obtained. What is required is the bulk heating of the loop plasma to 
a temperature of ^10® K as suggested by Chubb, Kreplin and Friedman 
(1968) to explain their early observations of hard x-rays of energies 20 keV 
and above. 

One type of mechanism suggested by Matzler et al (1978) involves the 
reversible adiabatic compression and heating of the plasma. There are, however, 
both observational and theoretical problems with this model (see Elcan 1980). 
Another type of mechanism is to accelerate electrons to high energies in the 
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reconnection region and to Inject these electrons as a high current beam Into 
the ambient flare plasma of low density. Plasma-beam Instabilities generate 
turbulence (Ion-cyclotron or Ion-acoustic* Kendell and Kennel 1971, also 
Duljveman, Hoyng and lonson 1981; Emslle 1981b) which enhances the electrical 
resistivity of the plasma and rapidly dissipates the electron beam depositing 
the beam energy as heat In the ambient plasma. The development of thermal 
models has been previously discussed (e.g.. Brown, Melrose and Spicer 1979; 

Smith and Mlliequlst 1979; Smith and Auer 1980; Emslle 1981a) and we shall not 
concern ourselves with such models here. 

The purpose of this series of papers is to analyze the general character- 
istics of the "non-thermal" models. All previous analyses of these models 

(e.g., Ramaty and Lingenfelter 1967; Lin and Hudson 1971; Brown 1972; Takakura 1972; 
Petrosian 1973; Kane 1974) have been limited to total flux calculations and have 

r 

either ignored the dispersion in phase space of the accelerated electrons or treated 
such in an approximate manner. The high spatial resolution capabilities of the 
HXIS experiment on the SMM and various on-going observational programs using the 
VLA and other high resolution microwave instruments will soon provide us with 
detailed hard x-ray images and inicrov/ave maps of flares during the impulsive 
phase. In order to utilize these observations we need to evaluate the distribu- 
tion of high energy electrons along the flare loop and calculate the resultant 
x-ray and microwave radiation as a function of position along the loop. 

In this paper we present the results of our analysis on the steady state 
distribution of energetic electrons within the flare plasma. We use the time 
independent Fokker-Planck equation to determine the steady state distribution 
in pitch angle, energy and height above the photosphere of the electrons as tliey 
spiral along the magnetic field lines and collide with the particles in the 
ambient plasma. In subsequent papers this distribution will be used to calculate 
the characteristics of the resultant x-ray and microwave radiation. In sll 
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we present the general features of this problem and the equations to be used 
In its solution. Some approximate analytic solutions are discussed in sill* 
and the results of the numerical analysis are presented in ilV. In IV 
we present a summary and our conclusions. 
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11. THE MODEL 

The model we are Investigating, which encompasses thick-target, thin-target and 
trap aspects of non-thermal hard x-ray nx)de1s. Is the following: High energy electrons 
with some Initial energy spectrum and pitch angle distribution are Injected at 
some height above the photosphere either Into a closed magnetic loop structure 
or into an open field structure (see Figure 1). These electrons, spiraling 
along the field lines, undergo collisions with the background particles (atoms. 

Ions and electrons) and these collisions, in addition to the variations In the 
magnetic field structure, diffuse the energetic electrons' distribution function 
in phase space. Note that wave-particle interactions, which may be important 
under certain circumstances, are not dealt with in this paper. In order to 
calculate the radiation from the electrons, we must first evaluate the distribu- 
tion in phase space of the electrons. 

t' 

Two characteristics of the background plasma simplify this calculation 
significantly. 

a) The time scale for energy loss by collisions is 

''con = - E0/.;,r2 ncloA (1) 

1 

1 2 2 2 

where E = y- 1 is the electron kinetic energy in units of m^c , p = 1-1/y , 

^All energies will be expressed in these units unless otherwise specified. 

p OA p 1 

4irr^ = 10" cm , and A is the minimum angle of deflection in the couloumb 

10 -3 

integral, typically InA - 20. At a density of 10 cm < 1 sec forel^trons 

!• 

with energy < 10 keV and is correspondingly longer for higher energy electrons. 
However, since the time L/pc for electrons to traverse a typical 
magnetic loop (L s 10^®cm) is less than a second, the more energetic electrons 
penetrate the high density regions and lose energy within less than a second 


/ 
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(cf. Petrosian 1973). It then follows that for the study of the Impulsive phase 
on time scales longer than this the electrons can be treated as If In steady 
state. 

b) For electrons with energies 10 to 1000 keV the gyro radius, 

Tg « 10 (100 gauss/B) cm, Is much, shorter than the characteristic length scales 
of the plasma, namely, the density or magnetic field scale heights, which are 

g 

Of the order of L s»10 cm. Contcquently, electrons ondifferent field lines 
are decoupled and the diffusion of the electrons perpendicular to the magnetic 
field can be ignored. Thus,, the electron distribution varies only along the 
magnetic field lines , which we take to be static. 

With these considerations the problem is reduced to a three dimensional 
time independent diffusion problem: One spatial dimension s along the field line 
and two momenta (parallel and perpendicular to the fie^ld lines respectively). 

It is more convenient to use as independent variables the cosine of the electron 

1 ■' ‘ ■'.» 

AAA A 

pitch angle, \i ^ p»b (p and B are unit vectors in the direction of the electron 
momentum and the downward magnetic field), the electron kinetic energy E and a 
dimensionless column depth t defined by 

dr = ds/X^, ^ 47rr^ n InA = 2 x 10"^\nr\lO^%"V)(lnA/20) , (2) 

v/here n is the number density of hydrogen and is a mean free path scale 
(the mean free path for electrons of energy E is X^Eb?-). Note that for InA = 20, 

T = 1 corresponds to a column depth of 5 x 10 cm which, in the quiet sun, 
is approximately where the temperature minimum is to be found. 

We use the Fokker-Planck method to solve for the diffusion in pitch angle and 
energy along a field line of given geometry and field strength. Let f(E,y,T)dEdudx 
be the number density of electrons v;ith cosine of pitch angle between p and 
p dp, and with kinetic energy between E and E + dE, lying between 
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column depths t and t + dt. The Fokker-Planck equation for the problem becomes 
(see e.g., Rosenbluth, McDonald and Judd 1957) 


8f 


2. dInB df 


JL 

dE 


>l> 




(3) 


where we have Ignored terms In the Fokker-Planck expansion which do not contain 
the large InA coefficient. The first term on the right hand side describes 

the mirroring of the particles owing to the convergence of the magnetic field. 

The last two terms account for the energy loss and pitch angle dispersion of the 
beam through collisions with the ambient particles. The coefficients preceding 
the energy dissipation term and preceding the pitch angle diffusion term 
depend upon the composition of the plasma. For a partially ionized plasma where 
X - n /n is the fraction of the hydrogen that Is ionized 


C = X .Ibxl InlsVE/a") 

h * ■ 2 InA • 



r 

1 


(4) 


where a is the find structure constant is the ionization potential 

for hydrogen in units of m„c ). For a fully ionized plasma x = 1 and the 

6 

coefficients and reduce to 

C, “ 1 . • <5) 

so that, for non-relativistic electrons, Cg ~ C^. 

Other than hydrogen the solar atmosphere contains only helium in signi^ 
ficant quantities. For a helium fraction, by number,' of 10%, the helium's 
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contribution to each of the coefficients and Is limited to 

s cJ/B, Cg® i Cg/Sg We shall Ignore the slight differences between the 
relative contributions of helium to and and shall ^\bsorb this 20 to 30 
percent contribution from helium Into the uncertainties In the values of InA 
and the density n. Therefore, we consider a pure hydrogen atmosphere. 

We have neglected the energy losses of the electrons due to the electric 
field needed to drive the reverse current, an assumption which holds well for 
high density and temperature plasmas and for small fluxes of accelerated elec- 
trons (tf. Knight and Sturrock 1977; Emslie 1980). The contribution of the 
energy loss due to the reverse current electric field can also be included in 
the above equation. However, this will complicate the solution considerably 
since the strength of the electric field at each level depends on the total flux 
of electrons at that depth F(t) = J Jp0cf(E,M,T)dEdp ,so that instead of a 
single differential equation we would have to consider a coupled integrb-* 
differential equation. This is beyond the scope of the present paper and 
will be considered in future work. 

In the next two sections we discuss some approximate analytic and some 
accurate numerical solutions of equation (3). 


V 
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III. SOME ANALYTIC SOLUTIONS 

In order to solve equation (3) we need to specify the characteristics of 
the plasma which determine the parameters C,, and dlnO/dt, and the boundary 
conditions. We assume that particles are Injected at the top of the loop 
(t » 0), such that the spectrum at t » 0, f^jlE.n) Is specified for 0 < p s 1. 
the other boundary condition being that f(E,M,T'H») » o. (The coefficients 
C| and Cg are approximately unity and have only a slight energy dependence.) 

The parameter dlnB/dt describes the variation of the magnetic field and deter- 
mines Its Influence on the electron beam according to the first term on the right 
hand side of equation (3). If the magnetic field and the density variation can 
be described by scale heights Hg and viz. n = n^, exp(s/H^); B = exp(s/Hg). 

then, in terms of our column depth t we find dlnB/dt = = H /tH„. Clearly 

o D n D 

when H^/Hg « 1 (I.e.. when the B field Is nearly uniform while the density 
varies rapidly) the magnetic field has an Insignificant effect and the electron 
distribution is modified by collisions alone. The resulting model Is a thick 
target model with the degree of beaming determined by the Injected spectrum. In the 
other extreme case, when H^^/Hg » 1 (I.e., low and nearly uniform plasma density 
and rapid field variation) the particles are reflected back and forth between 
mirror points of a closed loop configuration leading to a trap model with slow 
rate of precipitation, unless the density, n^. Is so high that the mirror points 
occur at T a 1. For an open field line configuration, particles are reflected 
no more than once and escape along the diverging field lines, giving rise to a 
situation as described by a thin target model (Datlowe and Lin 1973). However, 
it is usually found that the situation is more complicated owing to the variation 
of dlnB/dx along the field lines. 

In general, equation (3) cannot be solved analytically, and In the next 
section v;e shall describe a procedure for its numerical solution. Here we shall 
consider some special cases which allow us to derive exact analytic results. 
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The purpose of this is to check the accuracy of the numerical results and to 
derive some scaling laws which may simplify the presentation of those results. 

Case 1; If the plasma density is low and tha magnetic field scale 
height Hg is much smaller than the density height (so that 
dlnB/dx » X^/Hg » 1 throughout), then the last two terms on the right-hand 
side of equation (3) are negligible (C^bC 2 >> 0 ), and we have a situation corresponding 
to a collisionless trap model. It is then easy to show that the distribution function 
satisfies the differential equation 

8f/3t^ ® '-3f/aln(l-M^)» X ^ f dr “ In » ( 6 ) 

Jo 

giving a distribution function which is independent of energy and has the form 

f(T,u) =yn-(iv)e‘’)'*> . (7) 

where fjj(u) is the injected pi tch angle distribution'Lat t * o. This is simply 
a description of the adiabatic invariance of B/(l-p2). Note that for ari isotropic 
injected pitch angle distribution fjj(n) is independent of n and the pitch angle 
distribution remains isotropic throughout the atmosphere. For a gaussian injected 
distribution the distribution remains gaussian but with a width that increases with t. 

Case 2: If the magnetic field is approximately uniform such that 

Hn » H and X , then the first term on the right-hand side of equation (3) 

D n 0 

is negligible. However, even when dlnB/dr = 0, analytic solutions cannot 
be obtained except in the small pitch angle [(l-w) « 1 ] regime. Before presenting 
the solution for the general small pitch case, v/e first consider the follovn'ng two 
cases, which, though unphysical, expose the nature of the results for the physical 
case and which may be used for testing the numerical code (see Appendix). 

a) Eneray loss term with no pitch angle scattering (C^ = 1, Cg = 0, 
dlnB/dt =0). In this case it is useful to define 

4 . = f/B, dn 5 B^dE, n = E^/(E+1) 


( 8 ) 
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so that for a given (and constant) pitch angle y equation (3) and Its solutions are 

^aT'"'ln* n) » (9) 

where ♦J[n) Is related to the Injected spectrum f^CE) through definitions (8). 

2 

For non-relativl Stic energies n “ E « 3V4» and 

I 

Ne.iut) » f,(E/tT/iiE^) • 0 +t/ . (10) 

I 

(b) Pitch angle diffusion with no energy degradation (C^s*0, C^r » 
dlnB/dt = 0). In this case analytic solutions are possible only for small pitch 
angles. If we define 


«V2 5 1 - u « 1 and t ® C2t/3*tY^ • 


( 11 ) 


where a Is approximately equal to the pitch angle, then the differential 
equation (3) reduces to 


37 


1 3 f 3JF \ 

a 3a 3a^ 


which has the solution 


f(T,a) 


e”*^ ^A(w) J^(wa)do) , 


( 12 ) 


where Is the Bessel function of order zero and A(w) is determined by the boun- 
dary conditions. For example, if the injected pitch angle distribution is gaitssian, 
f (a) = (*~p-)e*‘“^'^“o v/ith a^ « 1 , then multiplication of both sides of equation (12) 

0 CL^ 0 

O 

by J^(u) a) and integration over all o for x = 0 gives 



substituting this Into equation (X2) wc obtain 






( 14 ) 


c) The general small pitch angle case with a uniform magnetic field 
(C^= 1. Cg® “ 0). In this case we have again oV2 “ (1 -m)« 1 so that 

with the help of definition (8) equation (3) Is reduced to 


4B**y^ 

3+y 




(15) 


which has the general solution (separable in t,i) and a). 


f(T,n»o) = 


/ 




j; 


B(\)dxj L‘c(n)] 

0 


tl)' 


J^((i)a)A(w)duj » 


where 


C(n) = exp 



( 16 ) 


(17) 


and B(a) and A (w) are determined from the boundary conditions. For an 

injected gaussian pitch angle distribution f(E,ii,T = 0) = f^(E)(2/a2)e““^^“o , 

where a^- « 1» from equation (8) we can obtain an expression for 
0 

(n) “ ff.iE(Ti)]/B(n)* Then, following procedures similar to those in deriving 
0 o 

equations (9) and (14), equation (16) generates 


f(E.n,r) = fj,(n+t) 


+ Inc 

o 


exp 


( o? + Inc 7 


(18) 


• c H ie(n+T)/C(n)l 
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In the non-relativistic limit c(n) « “ E** « p/V? and equation (18) becomes 


f(E,M.T) “ f^(c**E)c"**I2/(<»J+ lnc))e"*“*^^®‘o'*' , 


(19) 


C - 0 + t/E^). 

It Is normally more convenient to use the flux of particles, F(E,m,t), 

Instead of their number density, f(E,M,t), where the flux and number density are 
related by F(E,u,t) = cpf(E,n.T). The relationship between F(E,u,t) and the 
injected flux at t = 0, F^(E,»i) Is similar to the relationship between f(E,ii,t) 
and fjj(E,u) as given in equations (18) and (19), these relationships being 
obtained by replacing f and f^ by F and F^ and by changing p(n)/p(n+t) 
to lp(n)/p(n+T)l in equation (18) and c to c in equation (19). 

As is evident, the shapes of the pitch angle distributions at various t are 
identical provided the scaling parameter c is constant.' Similarly, for an injected 
power law flux spectrum F^(E) « E‘®, the” ratio of the energy spectra at various 
depths T to the injected spectrum have the same shape for constant values of c. 

In the next section we shall explore the extent to which this scaling law 
holds for large pitch angles and for non-uniform magnetic field configurations. 
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IV. NUMERICAL RESULTS 


We have solved equation (3) numerically using a modified version 
of the program developed by Walt et al (1968) for the related problem of 
auroral electrons (see also Walt 1967). For a prescribed Injection spectrum 
f(E,u) at T » 0, the program solves for f(E,u,T) for v > 0, utilizing 
either a reflecting or non-reflecting boundary condition at the top 
of the atmosphere (t» 0), thereby calculating the electron distribution 
function for a closed or open field line configuration respectively 
(see below). The details of the procedure are described In the 
appendix along with the results of our testing the numerical code against the 
analytic results described in the previous section. We now present some 
numerical results for the electron distribution function for the nine models 
described in Table 1. Wo use model 1 as our standard model and as a bat for 
comparison so that we may observe the effects of the parameters describing the 
ambient plasma, the magnetic field configuration and the Injected spectrum on 
the electron distribution function. 

Model 1 has an injected electron spectrum 

F(E,u>o*t»»o) " Fq(E,|i) ’ exp-|-“|”-| » (EO) 

with » 0.01 and fi 5 (a typical value from solar hard x-ray 
bursts), a fully ionized hydrogen atmosphere with a uniform 
vertical -magnetic field and an open, i.e, non-rofloctlng, boundary condition at 
the top. This means that electrons which retin'n to t « 0 and therefore have 
pitch angle cosines p < U allowed to exit from the top of the atmosphere 
and are lost from the calculation. In Figure 2 we show the evolution of the 
electron flux with depth as a function of pitch angle for this model. The 
evolution of the flux distribution is similar for different energies provided 
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the depth t is suitably scaled. For all except the highest energies, that is, 
at non-relativistic energies, the scaling expressed by the numerical results 
is the scaling t « E , Isee equation (19) with c " constant). 

Consequently, the curves in Figure 2 are labelled according to their value 
Of the parameter t/E . This scaling of t « E is to be expected as the electron 
mean free path is proportional to the square of the electron energy for non- 
relativistic energies. Near to the top of the atmosphere most of the electrons 
of a particular energy E were originally electrons injected at x - 0 with an 
energy slightly in excess of E and which have lost only a small fraction of 
their energy and which have been deflected only slightly away from their original 
pitch angle. As can be seen, the forward beam of electrons is still clearly 
defined and only slightly broader than the beam at injection. These curves 
also show that, even at small t/E , there is an, albeit negligible, component 

r* 

to the flux of electrons with u < 0. This effect is due to those few electrons 
which, despite the extremely low probability of being scattered through very 
large angles of the order of ir, have been scattered within the atmosphere suffi- 
ciently that they have returned towards the top of the atmosphere. In the 
process they have also lost a substantial amount of energy, indicating that 
these electrons originally possessed energies considerably higher than , 

E and were, therefore, injected in fewer number [see equation (20)). 

At successively deeper levels in the atmosphere, i.e., at larger values of 
t/E , the flux distribution as a function of pitch angle becomes increasingly 
more uniform and the forward beam broadens, reflecting the fact that at great 
depths the number of electrons with kinetic energy E to be found at any 
pitch angle comprises a very broad range of electrons injected at the top of the 
atmosphere with many different energies and with many different pitch angles. 
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At greater depths still, the total number of electrons at any energy falls (as 
Is to be expected), but the beam is still well defined^ 

The scaling described by equation (19), which was obtained from the 
small pitch angle approximation, appears to be valid through to much larger 
pitch angles. It breaks down at higher energies due to increasing relativistic 
effects (ci^ 1 +t/E } and the fact th'it our injected spectrum has an arbitra- 
rily imposed cutoff at an electron k'nvtia energy of 1 MeV. We also find that, 
in agreement with equations (18) and (19), the curves in Figure 2 may be approxi- 
mated by gaussians to a good .degree of accuracy even out to large pitch angles 
where the flux ratios have dropped to insignificant levels (a fraction of one 
percent of the forward, [n = 1], fluxes). The widths of these gaussians also 
agree with those expected from equation (19) as shov^n by the inset in Figure 2 
where the analytic gaussians (solid lines) are compared to the numerical 
values (points). * » 

In Figure 3 we show the energy spectrum at various depths for model 1 , 
with the ordinate being the integrated electron flux in the forward direction. 

At the top of the atmosphere the energy spectrum is the injected spectrum, E“®. 
However, as the beam moves through the atmosphere, the lower energy electrons 
lose energy more rapidly than the higher energy electrons and soon their. numbers 
become depleted, giving rise to M>e humped spectrum that forms throughout the 
atmosphere. 

Integration of equation (18) over pitch angle, dji = ada, yields, 
using the electron flux F(t,E) rather than the electron density f(T,E), 


F(t,E) « 



F„(n»-T) == F^(c*'E)/c’^^ 


(211 


which implies that, in the non-relativistic limit and for a power law injected elec- 

* p 

tron spectrum F (E) « E , the maxima of the above curves occur at E = t/6. As 
0 max 


17 


shown by the inset in Figure 3 this relationship holds well for numerical results. 

Furthermore, at non-relativistic energies and for a power law injected spec- 
trum. the flux in the forward direction normalized to that at the top, 

F(t.E,ii.1)/F„(E.n«l) » c * / (H«;*lnc) , (22) 

depends only on the parameter c. Figure 4 compares the values of this ratio with 
the values obtained from Ihe numerical calculations. For 30 keV and 70 keV elec- 
trons the agreement is excellent for the range of values of the parameter c perti- 
nent to solar flares. For 750 keV electrons we also compare the numerical 
results with the fully relativistic form of equation (18) and present this in 
the upper of the two curves for 750 keV. It is evident that the relativistic 
form indeed becomes necessary at high energies. Note that the 

* r 

analytic expression begins to deviate significantly (>20%) only at depths 
t/E > 0.3 where the electron flux is less than one percent of the injected 
flux (see Figure 2). Figures similar to Figures 3 and 4 may be drawn for the 
flux at any particular forward pitch angle but, given the good agreement between 
the results and the analysis as displayed in Figure 2, those figures obtained 
would clearly be very similar to Figures 3 and 4 which have been drawn using 
the integrated forward flux. 

Finally, if we integrate the above spectra over all depths, we obtain the 
spatially-integrated energy spectrum (which determines the radiation spectrum 
of the whole, unresolved loop). 


f^(F) -j F(T,E)dT =3 J Fp(n^T)dr = p?(n) j F^(E')dE' (23) 
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which* for E 


(E) 


E-«+'j2 . e -**2 


0 



( 24 ) 


Note that this result is not restricted to the small pitch angle regime. 
This can be demonstrated by integrating equation (3) first over m and then 
over T which gives 



= Fo(E) ^ 



+1 

F(E, M.x)dM 


(25) 


Integration of this gives equation (23) which v/as derived from the special case 
of p(6) “ 1 . . 

In Figure 5 we plot F(E) versus E for model 1 showing an agreement between 
our results and the above relationship which is especially good at low energies. 

“ 6+1 

This agreement breaks down at higher energies for two reasons: (i) F(E) E 
at extreme relativistic energies, and (2) the presence of an upper cutoff 
in the energy of the injected beam. 


Also shown in Figure 5 is the same function for models 5 and 9, both of 
which display a similarly good agreement with the analysis. Model 5 has a value 
of 6 = 3, and rnodeli' 9 is, parametrically, most distant from model 1 and from the 
approximations used in the analysis of model 1. 

In summary, the excellent agreement betv/een our numerical results and the 
above analysis of model 1, in addition to providing us with ample demonstration 
of the accuracy of the numerical results, indicates that the analytic results 
hold to a high degree of accuracy well outside the regimes of validity suggested 
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by the approximations used (that is, to large pitch angles) and that the 
analysis provides a good description of the evolution of the electron beam 
throughout the atmospiiere. 

Results for Other Models 

We now consider the results from the other models described in Table 1 
and compare them wi th those from model 1 . 

Model 2 has a strongly converging magnetic field throughout the atmosphere 
but retains the open boundary condition at t = 0. The effect of the converging 
magnetic field is to enhance the broadening of the pitch angle distribution with 
depth and to increase significantly the number of electrons which are moved to 
negative pitch angle cosines and return towards the top of the atmosphere. 
However, despite the presence of the magnetic field which converges by a factor 
of roughly thirty from t « 0 to t “ 1, the narrowness of the injected pitch 
angle spectrum is such that the downward traveling electron beam contains 
about two orders of magnitude more electrons than are removed from the beam 
by the magnetic mirroring so that the beam stays well defined throughout the 
atmosphere. 

Model 3 is as model 2 but investigates the effects of reflecting those 

electrons with n < 0 at t =0 back into the loop vnth u > 0. This 

mimics the effects of having the electrons in a closed and symmetric loop 

configuration (see Fig. 1). However, since the number of electrons exiting 

with n < 0 is relatively small, the effects of this change on the boundary 
condition at t - 0 is negligible. 

Model 4 is the same as model 1 except that it has an atmosphere 

19 -2 

which is fully ionized down to a column depth of 5 x 10 cm and 
then becomes fully neutral as rapidly as can be handled by the 
numerical code and assumes chromospheric densities throughout the 
remaining depths. Thus, model 4 introduces a transition region and has changing 
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values of the parameters and Cg [see equation (4)J. The effects, 
however, on the overall electron distribution, of having a partially ionized 
atmosphere, are small with the beam being only slightly affected helw the 
transition region (cf. Emslie 1978). 

Figure 3 shows the energy spectrum as a function of depth for model 1. 

A similar plot for models 2 through 4 shows that the energy spectrum at each 
depth is almost completely insensitive to the changes in the parameters introduced 
with these models and the humped spectrum is maintained with the same scaling 
En,ax® ■f- This is as expected because models 2 through 4 remain primarily 
collisional in their treatment of the beam despite the introduction of the 
converging magnetic field, etc., and therefore exhibit the same scaling as 
model 1. 

In general, we therefore conclude that for a narrow injected pitch angle 
distribution the change in the plasma parameters and the introduction of a 
reasonable magnetic field configuration have a negligible effect on the evolu- 
tion of the electron beam within the upper part of the target and tend to 
make the electron flux more uniform in the lower parts where the beam flux is 
less than 0.1% of the injected beam flux. 

Model 5 is identical to model 1 but has a harder electron 
energy spectrum with 6 == 3 (corresponding to the hardest observed x-ray 
spectra) in place of 6 = 5. Having more high energy electrons in the beam 
at the top of the atmosphere simply leads to having more electrons with medium 
energy at large pitch angles lower in the atmosphere and a larger number of 
electrons reflected back upwards with \\ < 0. Otherwise this model, like 
model 1, agrees well with the analytic expression given by equation (19) indjjsLated, 
for example, in Figure 5. 

Mociel 6 has a much broader injected gaussian pitch angle distribution; 

= 0.125. In Figure 6 we show the evolution of the electron flux distribution 
for this model where, comparing Figure 6 with Figure 2, we can see that the 
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effect of having a broader Injected beam is solely to give rise to a broader beam 

traveling downwards through the upper atmosphere. However, the shape of the 

flux distribution at greater depths becomes identical to the shape of the 

distribution function at the same depths in Figure 2. Remembering that the 

2 

curves in Figures 2 and 6 give the beam evolution with r scaling as E , 
we can see that for each energy there is a depth below which the shape of the 
distribution function in pitch angle is independent of the injected gaussian 
width. This is to be expected from equation (18) owing to the fact that Inc 
increases slowly with t until it eventually overwhelms o* even for the broad 
injected beam. Thus, below a certain depth the electron flux distribution 
relaxes to a smooth form and then slowly isotropizes. Only at exceptionally 
large depths where the electron flux is negligible does the distribution become 
sufficiently isotropic that the forward component is not easily discernable. 

The inset in Figure 6 is similar to that in Figure 2 apd shov/s the extent to 
which the small pitch angle approximation agrees with the numerical results 
for this model. Considering that the injected spectrum is broad, this agreement 
is better than might have been expected, the agreement holding well out to 
large pitch angles. 

Model 7 is the same as model 1 with the exception of now having an injected 
pitch angle spectrum which is uniform for positive y- The evolution of the 
pitch angle spectrum with depth is shown in Figure 7 and again we see that the 
shape of the spectrum at the two greatest depths is nearly identical to those 
in Figures 2 and 6 despite the extreme form of the injected spectrum. 

In models 8 and 9 the magnetic field is non-uniform (as in models 2 and 3) 
but v/ith an injected spectrum which is uniform as in model 7. In addition, model 9 
has a reflecting boundary at the top. Little needs to be remarked upon for models 
8 and 9 other than that, as expected, the effects of the vertical converging mag- 
netic field and the reflecting boundary condition are to give almost a symmetric 
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pitch angle distribution with f(v) “ f(*u) at each depth, the discrepancy being 

the leakage of electrons to lower energies by collisions. For comparison with 

Figure 3 the energy spectrum evolution with depth of i.iode1 9 Is shown In 

Figure 8. It can be seen that, even for model 9 which Is parametrically most 

distinct from model 1, the energy spectrum retains the same form with the peak 

In the spectrum moving to higher energies with Increasing depth In accordance 
2 

with T “ as shown In the Inset. 

maX 
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V. SUMMARY AND CONCLUSIONS 


1. If the Injected spectrum Is strongly concentrated In the direction 
of the magnetic field lines (small pitch angles), plausible variations in the 
parameters describing the plasma and the field configuration have significant 
effects only at large depths and large pitch angles where the beam density has 
fallen by several orders of magnitude. Otherwise, the analytic results of 
equations (18) and (19) give an acceptably accurate description of the evolu- 
tion of the beam. 

2. Even when the injected electrons are weakly beamed or not beamed at 
all (uniform pitch angle distribution), the scaling derived from the small 
pitch angle regime seems to be valid except at small optical depths (see 
Figures 6 and 7). 

3. At large depths the pitch angle distribution 'is independent of the 
injected distribution, evolving to a smooth and broad distribution. 

4. The energy spectrum evolves independently of the injected pitch angle ' 
distribution and agrees quite well with the analytic approximation given by 
equation (21). 

5. Most significantly we find that for most practical applications 
equations (18) to (24) may be used instead of having to resort to detailed 
numerical calculations. 

In future papers in this series we intend to utilize the results discussed 
and described here to calculate the radiation observed from a flaring loop. 

Using the complete electron flux distributiorf, we will be able to calculate the 
polarization and directivity of hard x-rays from the loop, the hard x-ray flu.x 
as a function of position along the loop and the details of the microwave flux 
including its dependence on the magnetic field structure defining the loop. 
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With results from current observational programs at the VLA (Marsh et a1 1980) 
and the high resolution data from HXIS (Hoyng et a1 1981). to mention but two 
of the many sources of high quality observational data, we hope to be able to 
use our results as a powerful diagnostic tool for the further understanding of 
the impulsive phase of flares. 

The model as it stands has several limitations which will receive atten- 
tion in the future pursuit of this work. The effects of reverse currents may, 
under certain physical conditions of the ambient plasma (see Emsliel980, 1981b), 
be of significance in determining the evolution of the flux of accelerated 
electrons in a flare limb and are intended to be included in the model to give 
it greater applicability for a broader range of flare parameters. Similarly, 
the effects of wave particle interactions on the energy spectrum at energies 
below that of maximum flux will also be investigated. ». 
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APPENDIX 

A. Description of the Numerical Code 

We consider a (horizontally uni form) I -dimensional plasma atmosphere with a 
density structure corresponding approximately to the low corona and chromosphere 
and Inject electrons at the top of the atmosphere, designated as t ** 0. 

The electrons are Injected according to a prescribed pitch angle distribu- 
tion with the electron's pitch angle being measured as the angle between the 
electron's Instantaneous velocity vector and the magnetic field line which Is 
in the downward direction along a flare loop (see Figure 1). Thus, downward 
moving electrons have positive pitch angle cosines; 0 < u .^1. These electrons 
then travel through the plasma atmosphere undergoing collisions with the cold 
plasma particles and being adiabatically scattered by the converging B field. 

They thus change their pitch angle and energy and are allowed to return to the 

f 

top of the atmosphere, then having negative pitch angle cosines. As the electrons 
can only diffuse from high to low energy and no process exists within our formu- 
lation for accelerating electrons to higher energies, the most convenient stepping 
parameter for the numerical analysis is the energy of the beam electron. The 
electron flux distribution is calculated over a three diminsional phase space 
(E,m,t) which is stratified into energy layers with each layer consisting of a 
grid in (u,t) space and being of constant value in energy space (see Figure A-1). 
The code generates the value of the electron flux at each grid point in a (p,t) 
plane of energy from the given input fluxes at t = 0 and the flux values at 

the overlying (m,t) plane of energy -j > E.. 

The input flux is specified at the top of the loop, t = 0, as a spectrum in 
energy and forward pitch angle, i.e., for the 0 < ji 1 snd E < 1 MeV half" 
plane. The flux values on the (w,t) grid con\..>pond1ng to the largest electron 
energy are just the input fluxes specified at the top of the loop for that 
energy (1 HeV). 


i 
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Once the flux values on one energy plane are calculated satisfactorily 
the code steps down one energy step to the next energy plane. The fluxes from the 
preceding plane just calculated are scaled for the new energy and, along with the 
specified Input flux at that energy, a zeroth order (h,t) grid of fluxes is gener- 
ated. The code then Iterates over this energy plane until a satisfactory conver- 
gence of the flux values to those values satisfying equation (3) Is found. A 
final n^^ order t.,T)gr1d of electron fluxes Is obtained which Is stored as part 
of the results and Is used to generate the zeroth order gr^j for the next 
energy level down. Iteration on each (n.T)grld Is obtained by sweeping down from 
the top of the atmosphere (t*" 0) to the bottom and then returning to t - 0. 

The process Is repeated until the results after each pair of sweeps have converged 
to v/lthin a prescribed error. 

Boundary conditions can be imposed at the top to allow for electrons 
exiting with u < 0 to either be lost from the region of the calculations 
(simulating an open field configuration) or to be reflected back with 
positive Vi's to simulate a closed field configuration symmetric at t = 0 
(see Figure 1). Electrons with u > 0 at the bottom of the atmosphere are 
lost as are electrons which have their energy degraded to less than the minimum 
energy, which we have set at 10 keV. 

^ • T ests of the Numerical Results 

In order to assess the accuracy of the code, v/e compare the numerical 
results with those of the analytic results described In sill. Note that 
some of the cases considered here have no physical meaning and servo solely as 
examples for which the analysis is tractable. 

Case 1 ; If we set the coefficients C-j and C 2 in equation (3) equal to 
zero and inject electrons v/ith a pitch angle distribution f^Cvi) = -k1n(l- vi^O* 
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then, according to equation (7) we expect the pitch angle distribution to vary 
with depth In the form f(T*ij) « k{ 7 - In(l-y^)) , that is, curves of f(T,u) 
versus -ln(l-M^) at constant 7 , or depth, will be straight lines which are 
parallel for different depths and vertically positioned according to the value 
of Ict at each depth. In Figure A-2 we have plotted our numerical evaluation 
of f(T,u) against corresponding values of -ln(l-y2) for suveral depths. As 
the electron's pitch angle Is being affected adiabatically by the magnetic field, 
there is no energy dependence or energy loss In the process by which the electron 
beam evolves and, therefore, Figure A-2 is identical for all electron energies. 

As can be seen from Figure A-2, the numerically evaluated f(T,u) behaves 
exactly as expected throughout the atmosphere except at large depths where the 
difficulty of prescribing the injecteti flux In finite sized pitch angle bins 
for values of u close to unity (where -ln(l-n^) becomes very large] becomes 
evident. ^ 

Case 2a ; If we now set = 1 and C2 “ dlnB/dr = 0, we mimic a situa- 
tion where the electrons lose energy by collisions within the atmosphere but do 
not change their pitch angle. If we inject a spectrum = f^(E)/3(E) - 

[E^/(E + 1)]”^, i.e., a pov/er law in n (recall that n = E^(E + 1)], then, 
according to equation (10), we expect the electron distribution to be described 
at each depth by (|>(E,u,t) - (t/u + n)~*. In Figure A-3 we have plotted 
{<|i(E,p,t)“^^* - 4>(E,p, 0)'^^^)for three values of p, versus depth t and note 
that this should be independent of E and proportional to t/p. As is evident 
the numerical results agree perfectly v/ith the expected analytic relation. 

Case 2b : Finally, to check the accuracy of the treatment of diffusion 

in pitch angle we set C^ - dlnB/dr = 0, C^ = (3 +y)/4 and inject into the 
atmosphere a beam with gaussian pitch angle distribution. Equation (14) 
describes the evolution of such a gaussian beam but is an approximation valid 
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at small pitch angles. Plotted against gausslan distributions at each 
depth, as given by equations (11) and (14), would be straight lines of slope 
■•l/(aj + 4t), and these are shown as solid lines In Figure A-4 (for o* ® 0.01). 
The points and broken lines are the results from our numerical analysis at 
several depths t. As can be seen, for small pitch angles the numerical 
results agree with the expected result to a high accuracy. Furthermore, as 
discussed In SIV, we note that the small pitch angle approximation Is 
valid to larger angles than expected, especially at higher depths. At every 
depth the approximate analytic expression deviates from the numerical result 
only when the fluxes have dropped to less than one percent of the fluxes in the 
forward direction. 
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FIGURE CAPTIONS 

Figure The two magnetic field structures that we consider. One is an 
open field structure with electrons injected (at t = 0) at some height 
above the photosphere with reflected electrons being allowed to leave 
freely once they reach t » 0 with u < 0 . The other structure is the 
closed symmetric field structure where the electrons are injected at the top 
of the loop (t = 0) down towards the photosphere symmetrically in both 
directions. The electrons which return to t = 0 with u < 0 are then 
reflected back into the volume of positive t to simulate the electrons 
crossings from one leg to the other. 

Figure 2. The form of the electron flux pitch angle distribution at several depths 
in the atmosphere for a narrow injected beam (model 1) as obtained from 
the numerical results. Electrons with a pitch angle of zero degrees are 
streaming along the magnetic field linos downwards towards the photosphere; 
those with a pitch angle in excess of 90 degrees are moving upwards. 

The depths are scaled with the square of the electron energy as explained in 
the text. The inset shov/s a comparison of the numerically derived values 
(points) to those expected from the analytic expressions of equation (18) 
(solid lines). 

Figure 3. The form of the electron flux energy spectrum at several depths in the 

atmosphere for model 1. At injection the beam has a flux spectrum of the 

form but low energy electrons lose energy by collisions faster than 

high energy electrons and are rapidly depleted in number v^fithin the atmosphere. 

The spectrum then develops a hump which moves to higher energies at greater 

2 

depths according to E,„.,^ = t/6 as shown in the inset. 

lllCl ^ 
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Figure 4. A comparison of the numerical results with the analytic expression 

of equation (22) at three energies. The upper diagram which is for energies of 
30 and 70 keV uses the non-relatlvlstlc expressions for c which should 
be a good approximation at these energies. The lov/er diagram Is for 
an electron energy of 750 keV and shows a comparison of the numerical 
results with two forms of equation (22), the lower one using the non- 
relativistic expression for c and the upper one using he relativistic 
expression for c* 

Figure 5. The energy spectrum of the spatially integrated electron flux for three 
of our models: model 1 (our standard model); model 5 ( which is the same as 
model 1 except that it has a harder injected electron energy spectrum); and 
model 9 (which is the most different from model 1 [see Table Ij). According 
to the analysis the total electron flux energy spectrum should have a partic- 
ularly simple form as given by equation (24) and Ishown by the solid lines. 

The points represent the numerical results for the three models. 

Figure 6. Same as Figure 2 i but for model 6. 

Figure 7. The same as Figure 6 but for a beam which is isotropic in pitch 
angle at injection (model 7). Though the injected beam is not 
a gaussian, the form of the pitch angle spectrum at the two greatest 
depths is the same as those for model's 1 and 6 : 

Figure 8. The electron energy spectrum at several depths in the atmosphere 
for model 9. As can be seen from Table 1, model 9 is most distinct from 
model 1 in terms of its parameters, yet, as a comparison with Figure 3 will 
show, model 9 retains the same basic form for the energy spectrum as model 1. 
This is to be expected as, for both models 1 and 9, the electron beam is 
primarily coll isionally modified. 
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Figure A-1. A schematic respresentatlon of the three dimensional phase space 
of the electron number distribution. The phase space is stratified Into 
(horizontal) energy planes, with each plane being a grid In pitch angle 
and depth within the atmosphere over which the electron number distribution 
Is evaluated. The electron flux Is specified over the positive pitch 
angle half of the t ® 0 plane surface (shown shaded), the lower boundary 
being put at t = « for which the number distribution is taken to be 
Identically zero. The other half of the t = 0 surface (unshaded) corres- 
ponds to electrons with y < 0 at the top of the atmosphere, and these 

electrons can be either allowed to leave this region of phase space and 
lost from the calculation (open field configuration of Figure 1) or be 

reflected back Into the region with y > 0 and added to the already specified 

Injected flux (closed loop in Figure 1). 

Figure A-2. The results of testing that part of the code handling the inter- 
action of the electron beam with the converging magnetic field. The solid 
lines show f(y,T) as given by equation (7) with fjj(y) = -kln(l-y^); 
the points represent the numerical results. The right hand side of the 
diagram is for downv/ard travelling electrons, the left hand side for 
reflected electrons travelling upward. The depth parameter, 7 , is related 
to the convergence of the magnetic field in the form 7 = ln{B(x)/B(0)). 

Figure A-3. The results of testing the part of the code describing the energy 
degradation only for electrons with pitch angle y = 1.000, y = 0.9563 
and y = 0.788. The solid lines represent equation (9) for an injected 
electron spectrum ({.^(E.y.r) = (E /E+1) ; the dots represent the 

numerical results. 


Figure A-4. The results of testing that part of the code dealing with the 
diffusion of the electron's pitch angle by collisions for an injected 
gau&sian pitch angle distribution. Note the excellent agreement between 
the numerical results (points) and the analytic expressions (solid lines) 
for small pitch angles. 
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